* ANALYSIS AND DESIGN OF MODEL FOLLOWING CONTROL SYSTEMS BY STATE SPACE TECHNIQUES * 
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control system and of deciding when the plant 
cjn follow the model without error is consid- 
ered. Specifically, a testing procedure is 
given which determines when the output of a 
linear plant can he made to follow a model by 
placing feedback around the plant. It is 
assumed that model and plant are described by a 
known set of state equations. The tests are 
developed for two configurations of model fol- 
lowing. In the first of these, known as 
implicit model following or matching dynamics , 
the model does not actually become a part of 
the total system but enters only into the selec- 
tion of the feedback law, whereas in the second 
or real model following configuration, simula- 
tion of the model is required. It is shown 
that in the absence of disturbances the condi- 
tions for perfect following are essentially the 
same for both configurations. For the case of 
implicit model following, the control law which 
achieves perfect following is also calculated 
and, in general, will require both finite and 
singular controls of the delta function type if 
the model state space is of smaller dimension 
than the plant state space. The results 
obtained suggest a rational guide for deciding 
what type of model following is most appropri- 
ate for a given problem. That is, real model 


methods of deriving an adequate model but 
assumes at the outset that one is given as part 
of the problem. 

Recently Ellert and Merriam [l] and Tyler 
[2] used quadratic optimal control theory to 
synthesize model following control systems. 

Their technique, unlike those based on classi- 
cal procedures, is applicable to arbitrary 
multivariable systems and always yields a feed- 
back configuration which minimizes a quadratic 
function of the error between the plant output 
and the model. 

Although quadratic optimal control yields 
a general procedure for synthesizing model fol- 
lowing control systems , experience has brought 
to light some additional design problems which 
can be handled more effectively by other meth- 
ods. An example of such a problem is that of 
deciding when the closed -loop plant can follow 
the model perfectly. That is , for a particular 
combination of plant and model, one may find 
that the closed -loop plant designed by the meth- 
ods of optimal control follows the model with 
unacceptably large errors which cannot be 
reduced below some limiting value merely by 
manipulating the weight matrices in the cost 
function. In other words, there may not be 
enough freedom of choice in the feedback matrix 
to match the plant to the model if the model 


following, at the expense of greater complexity, dynamics differ greatly from the open-loop 


offers the best performance if random distur- 
bances occur within the plant but achieves no 
better performance than the simpler implicit 
model following in the absence of disturbances 
or uncertainties. Two examples, one of which 
is based on the lateral equation of motion of 
an aircraft , are given to illustrate the theory. 


Introduction 


plant dynamics. In that case, the feedback 
matrix calculated via optimal control still 
yields a least squares match between model and 
plant response during the control period but 
gives no prior indication of matching accuracy 
which must be determined separately either by 
actually checking the response of the closed 
loop system or by evaluating the minimum cost. 

A simple test of whether a plant and model can 
be matched when no restrictions are placed on 
the magnitude of the control signals would 

The design of a model follower control sys- therefore be a useful aid in designing a model 
tern consists of choosing a feedback law so that following control system, 
the output variables of the plant will faith- 
fully follow the output variables of a model. 

For instance, the plant could be represented by 
the linearized equations of motion of an air- 
craft and the model by the equations correspond- in the performance index" or implicit model, 
ing to an aircraft with ideal response charac- The relationship between these two design 
teristics as determined perhaps from simulator approaches is not entirely clarified in the lit- 
studies . This pape r is not concerned with thfc ePature , although Tyler [2] has shed some light 
*0ctober 18 , 1967 on this question. 


A further difficulty occurs when faced 
with the choice between a design based on a 
"model in the system" and the so-called "model 


' The main subject of this paper is to 
develop simple algebraic tests, applicable to 
both types of model following designs, for 
checking whether a plant can be matched per- 
fectly to a model. The tests are of such a 
nature that they also reveal the necessity of 
an unbounded control law if required for per- 
fect matching. For implicit model following, 
if the control law which achieves matching is 
known to exist, it will also be derived. 


Implicit Model Following 

In implicit model following one attempts 
to modify the output dynamics of the plant by 
means of feedback so as to approximate the 
dynamics of a given model. Ifethematically, 
implicit model following is defined in the fol- 
lowing way. Let the multivariable plant be 
described by the equations 

x = Fx + Bu (l) 

y * Hx (2) 

where x is an n -dimensional state vector, u 
is an m-dimensional control vector, and y is 
an Z -dimensional output vector. The matrices 
F, B, and H do not depend on time and have 
dimensions nXn, nXm, and ZXn, respectively. 
Also, it is assumed that n ^ m and n ^ Z. If 
the model is described by the equation 

z = Lz L = Z X J constant matrix ( 3 ) 

where z denotes the Z -dimensional state vec- 
tor of the model, then we ask that the output 
y(t) satisfy 

^ = Ly(t) 00 

as closely as possible. One method of achieving 
this objective is to use optimal control theory 
to minimize the following quadratic loss func- 
tion [3] : 

p 00 

I = / [(y - Ly)’Q(y - Ly) + u*Ru]dt 

^ o 

where Q is a positive semidefinite matrix and 
R is a positive definite matrix- This formu- 
lation of model following does not introduce 
the state variables of the model directly, 
since y and y appearing in the loss function 
can be expressed as a function of x and u. 
Hence, the names "implicit model," "model in 
the performance index" [2], and "matching 
dynamics" [4] all have been used to describe 
this method. Although optimal control theory 
generates a feedback matrix for arbitrary L, 
the error (y - Ly) obtained from the 


closed-loop system may be large even if the 
norm of R is chosen very small. Furthermore, 
one can show that assumptions such as control- 
lability of the plant or observability of the 
output also are neither necessary nor sufficient 
to guarantee that the error will be small. 

We now address ourselves to the question, 
"Under what conditions is it possible to sat- 
isfy equation (4) exactly?" Using equations (l) 
and (2) and requiring that equation (4) is a 
strict equality permits us to write 

HBu = (LH - HF)x (5) 

If this equation is to hold for all x by 
proper choice of u, then the range of HB 
must contain the range of (LH-HF) written as 

R(HB) Z> ft(LH - HF) (6) 

First, equation (5) is formally solved for u 
by taking the pseudoinverse of HB [ 5 ] 

u = (HB) f (LH - HF)x (7) 

Then, when u is eliminated from equations (5) 
and (7)> the condition for zero error becomes 

((HB)(HB) t - I)(LH - HF)x =0 all x (8) 

To justify the use of the pseudoinverse, one 
must show that if equation (8) is true for all 
x (i.e., ((HB)(HB)t - l)(LH - HF) is the zero 
transformation), then relation (6) is a neces- 
sary consequence. A property of the pseudoin- 
verse which permits this conclusion is that 
(HB)(HB)+ is an orthogonal projection operator 
on R(HB). Let z be any vector in ft(LH - HF) 
and write z as the sum 

z = z 0 + z ± 

where 

Zq e R(EB) and Zj _ e [R(HB)] X 

Since ((HB)(HB) t - l) is also an orthogonal pro- 
jection which projects every z e ft(LH - HF) on 
[R(HB)] and since by assumption equation (8) 
is zero for every x, it follows that 

z = z Q 

and therefore ft(LH - HF)cZR(HB). Thus, we 
conclude that choosing 

u = (HB) t (hH - HF)x (9) 

when 

((KB)(HB) + - I)(LH - HF) = 0 (10) 


'guarantees that y = Ly, or, equivalently, that 
the output dynamics of the closed-loop system 
will match the desired output dynamics. Fur- 
thermore, the boundedness of the pseudoinverse 
implies that the feedback law (HB)^(LH - HF) is 
bounded. Therefore, if the condition for zero 
error (eq. (10)) is satisfied and the model is 
stable, the controls that achieve a perfect 
match are always bounded. 

When equation (10) is not satisfied, it 
may still be possible to achieve zero error by 
enlarging the class of controls to include 
delta functions. As the next step the control 
law and the test for perfect following derived 
above are extended to the case of unbounded 
'controls. One begins by writing every control 
as a sum of ordinary and delta functions: 

u(t,r) = ui(t) + ug(r)5(t - r) (ll) 
me [3t(HB)] x 
ug€ 3l(BB) 

where !Jt denotes null space, J. denotes the 
perpendicular complement, r is a running vari- 
able, and t is current time assumed to be 
fixed. Thus the delta function occurs at time 
t. It is necessary to restrict ug to the 
null space of HB since otherwise the left 
side of equation (5) would contain a delta 
function of strength HBug while the right 
side does not. Hence, perfect matching would 
be absent at the moment the impulse occurs . 
From this remark it also follows that delta 
function controls are only helpful if the rank 
of HB is less than maximal. Adding the pro- 
posed delta function control at time t + can 
then be shown to modify the derivative of the 
output, y(t + ) , as follows 

y(t + ) = HFx + HBui + HFBug (12) 

Since the right side of equation (12) must be 
equal to Ly if zero error is to be achieved, 
one obtains 

HBui + HFBug = (LH - HF)x (13) 

By defining fl = ui + tog and using the pseudo- 
inverse to construct the necessary projections, 
one can write 

ux = (HB) f (HB)u , ug = (1-(HB) + HB)u 

( 1*0 

Then, upon substituting equations (lA) into 
equation (13), it is possible to solve explic- 
itly for ft: 

u = M + (LH - HF)x (15) 


where 

M = HB + HFB(l - (HB) + HB) 

Finally, the condition for zero error can now 
be derived by replacing Ux and ug in equa- 
tion (13) with the relationship for these 
quantities obtained from equations ( lA) and 

(15): 

(MM + - I)(LH - HF) * 0 (16) 

If condition (l6) is satisfied, then equa- 
tion (15) essentially gives the control law 
which achieves zero error, except for the 
implementation of the delta function control 
of that component of u which lies along 
9t(HB). Assuming for the moment that it is pos- 
sible to generate the required delta function, 
we want to demonstrate that from t + onward, 
equality of equation (A) can be maintained. In 
general, equation (A) or (5) will not hold at 
time t since the effect of the delta func- 
tion is not felt until time t + . At that 
moment a step change occurs in y in such a 
way that equation (A) is satisfied. Perfect 
matching is, therefore, as sirred for at least a 
time interval that is short in comparison to 
the fastest time constant of the system. As 
soon as the difference between y and Ly 
exceeds some small threshold, where the value 
of the threshold may be chosen arbitrarily 
small, another delta function whose weight is 
chosen according to equations (lA) and (15) is 
applied. The second delta function restores 
the equality of equation (A). Clearly, perfect 
matching can thus be maintained indefinitely by 
continuing to apply a delta function whenever 
the threshold value is exceeded. We also note 
that the smaller the threshold value is chosen, 
the closer will be the spacing of the delta 
functions, but also the smaller will be their 
strength. 

The problem of implementing a closed -loop 
control law which generates the required delta 
functions is discussed in the appendix. It is 
shown there that an approximate synthesis of 
such a control law is obtained by multiplying 
ug by a large positive gain constant K and 
that the approximation to the ideal delta 
function control law improves in proportion to 
the magnitude of K. 

If equation (16) is not satisfied and the 
rank of M is not yet maximal, it may still be 
possible to achieve zero error by including 
various derivatives of delta functions in addi- 
tion to the previously used controls. Consider, 
for example, the addition of first derivatives 
of delta functions. The control u is then 
written as the direct sum of three controls: 



t) (17) 


( 20 ) 


u = ui + Ug8(t - r) + UgiS 1 (t - t) (17) 
Ui€ [K(HB)] X 
u 5 e gi(HB) f\ [9l(HEB)] x 
u &1 e 3l(HB)f\ 3l(HFB) 

This decomposition assures that delta functions 
and their derivatives do not appear in the 
expression for y which, for the choice of 
controls given by equation ( 17) becomes 


y = HFx + HBui + HFBug + 


HF^Bi 


V 


At this point, the right side of the above equa- 
tion is equated to LHx and the condition 
describing when the resulting equation has a 
solution u for all x is found. As before, 
a numerical test for perfect matching similar 
to equation (l6) can be developed by defining 
u=ui+us+Ug 1 and then writing the compo- 
nents of u as° orthogonal projections of u 
on the appropriate subspaces: 


Ui = Piu , U5 = Pgu , Ug! 


= v* 


For general vector controls which contain 
finite, delta function, and derivative of delta 
function components it does not appear to be 
possible to construct the projections explic- 
itly in terms of pseudoinverses, although numer- 
ical procedures for performing such construc- 
tions are well known. The condition for 
perfect matching and the corresponding control 
law for this case is still given by equa- 
tions (l6) and (15), respectively, but M must 
now be replaced by 


M = HB + HFBPg + HF 2 BP, 


S 1 


(18) 


Real Model Following 


In implicit model following design the 
model entered only in the selection of the 
feedback matrix, and no real time error mea- 
surement between model output and plant output 
was necessary. In real model following the 
model, although itself uncontrollable, becomes 
part of the system in that the model states are 
compared with the output of the plant. The 
comparison actually takes place in the perfor- 
mance measure as follows : 


p O0 

I -J [(y “ z)*Q(y - z) + u'Ru]dt 


(19) 


where 


p 00 

I = / (w*Qw + u*Ru)dt 
o 

p °i, q.[ h,qh - h ' q i 

L 0 lJ L -qh q J 


• _ x 
w = -r- 

z 


It has been shown that the feedback matrix com- 
puted by this method depends only on the F, Q, 
and R and not on the model parameters L, 
whereas the feedforward matrix depends on both 
model and plant parameters [2], 

The derivation of conditions for perfect 
following in this case uses the fact that if 
all orders of time derivatives of the error are 
zero at any time t, then the error will be 
zero for all time. Beginning with the zeroth 
derivative one obtains the obvious fact that 


z(t) = Hx(t) 


( 21 ) 


If the first derivative of the error is to be 
zero, one finds after using equations (l), (2), 
(3)> and (2l) that 


H(Fx + Bu + BAz) = H(Fx + Bu + BAHx) = LHx 


If optimal control theory is to be used to com- 
pute the feedback and feedforward gain matrices, 
one augments the state space of the plant with 
the model states and then minimizes 


( 22 ) 


Here A is assumed to He an arbitrary feedfor- 
ward matrix. Upon solving the last two members 
of equation (22) for u, one obtains 

u = (HB) f (LH - HF - HBAH)x (23) 

This control law achieves zero error in the 
first derivative if 

((HB)(HB) t - l)(LH - HF) = 0 (24) 

which is obtained by substituting u of equa- 
tion (23) into equation (22). But equa- 
tion (24), if true, holds for arbitrary t; 
therefore, all higher order derivatives of the 
error will also be zero. Thus, the condition 
for perfect following for this case is the same 
as for implicit model following; moreover, 
direct use of the model states z(t) through 
the feedforward loop has no effect on this 
condition. 

If condition (24) is not satisfied, then 
one can consider sums of ordinary and delta 
function controls as a means of achieving zero 
error, just as for the case of implicit model 
following. Clearly, the arguments presented 
there carry over to this case. 


Evaluation of Real and Implicit 
Model Following 

It has been shown that the conditions for 
perfect following and the control law that 
achieves perfect following are identical for 
both real and implicit model following. Thus, 
assuming that perfect following is possible 
with either bounded or unbounded controls and 
that unknown disturbances are absent, there is 
no essential advantage of one design over the 
other. The key issue in deciding between a 
real model following design (with its addi- 
tional hardware requirements) and the simpler 
implicit model following is whether or not the 


designs , one calculated by the theory developed 
in this paper, and the other two by the methods 
of quadratic optimal control, are compared in 
this example, and the advantages of each are 
pointed out. The second example illustrates 
the theory when an unbounded control law is 
required for perfect matching. All computations 
were performed with the automatic synthesis 
program of Kalman and Englar [4]. 

Example 1 

The numerical values for the model and 
plant parameters used here correspond with one 
of Tyler’s examples [2]. 


requirements of the problem dictate that a par- 







ticular phase trajectory of the model be fol- 
lowed in the presence of unknown disturbances 
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in the plant. Implicit model following is not 


0 

-2.93 


-4.75 

-0.78 

capable of following a phase trajectory of the 
model where disturbances are present since no 
real time error measurement between model and 

F = 

0.086 

0 


-0.11 

-1.0 

plant states takes place; the model following 


0 

-0.042 

2.59 

-0.32 

is open loop as it were. But, if the model 







serves merely to characterize the desired 
dynamic properties of the plant , in other words 


~o 

0 




model and plant should have similar responses 
when starting at the same initial states with 

B = 

0 

-3-91 




no disturbances present, the implicit model 
following would be sufficient. 


0.035 

0 






-2.53 

0.31 




The maintenance of alinement between plant 







and model in the presence of uncertainties, be 
they unknown parameters or random disturbances 
necessitates the use of a real model in the 

H = I 

"o 

l 


0 

0 “ 

system. With a model in the system, errors 
arising between model and plant states due to 
uncertainties can be measured and corrected 

L = 

0 

-1 


-73-14 

3.18 

continuously. Thus, the principal advantage of 
having a model in the system is not that it 


0.086 
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-0.11 

-1.0 

always achieves better following, but that it 
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0.086 
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desensitizes the following to unknown 


1 — 




disturbances . 



cpj (bank angle) 



In the case of real model following, the 
control law given by equation (23) cannot be 
used by itself since it does not include the 
states of the model. That is, this control law 
fails to take advantage of the possibility, 
unique to real model following, of realining 
the plant and the model states if disturbances 
cause them to drift apart. Here the techniques 
of optimal control would seem most appropriate 
for computing the control law. 


Examples 


State vector 


Control vector 



(bank rate) 

(sideslip angle) 

(yaw rate) 

(rudder deflection) 
(aileron deflection) 


The test of perfect following, equation ( 8 ), 
applied to this example gives the following 
result: 


In this section two examples are presented. 
The first, which is discussed in some detail, 
represents the linearized lateral equations of 
motion of an aircraft. Three model following 


((HB)(HB) t 

- I)(LH - HF) 

= 


_ 0 

0 

0 

0 

-1.3X10 -7 

-4.3X10" 6 

-1.1XL0" 5 

-3 . 4X10' 6 

-1.2X10" 4 

-3.9X10" 3 

-1.3X1C 2 

-3.0X10' 3 

-1.6X10" 6 

-5 . 4xi0 -5 

-1.8X10 -4 

-4. 0X10 ' 5 


(25) 


Thus perfect following is not possible because 
the right side of the above equation is not the 
•zero transformation. Since HB has maximum 
rank, it also follows from earlier work that 
delta function controls cannot improve this 
situation. Nevertheless, because most entries 
in the matrix of equation (25) are small in 
comparison with entries in the system and model 
matrices, it is interesting to compare, the 
performance of the simple model following con- 
trol law of equation (7); that is, 

u = (HB) f (LH - HF)x = 


-3. 4X10 “ 3 
0 


-0.11 -0.37 

-0.49 17-5 


-0.084 


x 


- 1.01 


(26) 


with those calculated by optimal control for 
both the implicit and real model following per- 
formance indices. The Q, R feedback and 
feedforward matrices used in the calculation 
for the implicit and real model following 
designs are given below: 

Implicit Model Following 
Diag Q = [0,6, 0,6] , Diag R = [1,1] 


Feedback matrix = 

0.0034 0.111 0.371 0.0356 

0. 0.494 -17.5 0.6l4_ 

Real Model Following 

Diag Q = [10,10,10,10] , Diag R = [l,l] 
Feedback matrix = 

-0.074 -0.094 2.34 -3.23 

-3.15 -2.73 0.835 0.261 


Feedforward matrix = 

~0.031 -0.246 2.109 -4.62" 

-3.02 -3.16 1 6.6 -2.3 _ 

Figures 1 and 2 compare the transient 
responses of the three different control laws 
for two initial conditions corresponding to an 
initial bank angle and an initial bank rate, 
but with the model and plant states alined at 
the start. Because the response of the 
implicit model following law calculated with 
optimal control was generally not much differ- 
ent from the response obtained by using the 
control law equation (26), it is not drawn in 
all the figures in order to reduce crowding of 
the curves. Also, those state variable time 
histories that were omitted were found to be as 
well matched as cp in figure l(a). It can be 
seen in figures 1 and 2 that at least during 
the first 5 seconds , the performance of the 
control law given by equation (26) compares 
favorably with both real and implicit model 
following designed via optimal control. When 
model and plant are not too dissimilar, as 
shown in this case by the results of the per- 
fect following test, one can expect this con- 
trol law to work quite well; but for greatly 
mismatched model and plant, again as determined 
by the perfect following test, no assurance of 
satisfactory operation can be given. 

Figures 3(a) and 3(b) demonstrate when it 
is advantageous to use real model following. 
Here a disturbance in the plant is assumed to 
have caused a sudden misalinement between the 
model and plant bank angle variables. Under 
this condition, the real time error measurement 
between model and plant, which is only possible 
with real model following, facilitates the 
eventual realinement of corresponding state 
variables . Thus , the model serves as a memory 
of a particular trajectory in the presence of 
disturbances . 

Example 2 

The open loop plant equations of this 
example are again fourth order but the model 
equations are now of second order. 



Conclusion 
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Two basic approaches , each having its par- 
ticular advantages and disadvantages, exist for 
designing a model follower control system. In 
the implicit model following method the model 
enters only into the selection of the feedback 
law placed around the plant but does not become 
physically part of the total system. Thus, 
feedback is used chiefly to modify the dynamics 
of the plant so that its output behavior coin- 
cides with that of the model. This type of 
following therefore operates open loop with 
respect to the model since during the control 
interval no real-time comparison of model 
states and plant output takes place. The main 
advantage of this method is simplicity and low 
cost of implementation because the model in the 
system need not be simulated. 


It can be shown that the test for perfect fol- 
lowing with finite controls, equation (10), 
applied to this example fails; therefore, it is 
necessary to use the more general test given by 
equation (l6), which considers sums of finite 
and delta function controls. This latter test 
shows that perfect following is indeed possible 
with a control law containing both finite and 
delta function controls. 


After the required calculations are per- 
formed, the two parts of the control law are 
found to be 


ui = 


0 1.465 2.875 0\ 

.0 1.465 2.875 0/ 


If design specifications require that the 
model follower control system be able to follow 
a specific phase trajectory of the model start- 
ing at a given initial state while the plant is 
subject to unknown disturbances or parameter 
changes, then real model following is the appro- 
priate choice. Here the continuous measurement 
of error between model states and plant output 
offers the additional freedom of using this 
error, appropriately weighted, as a means of 
alining the model and the plant. However, 
because the conditions for perfect following 
are identical for both real and implicit model 
following, this additional freedom does not 
contribute to improved matching of the dynamics 
of model and plant in comparison with implicit 
model following. 


-0.086 -2.0 


-1.89 


u 6 = l x 

0.086 2.0 I.89 -1/ 


As shown in the appendix, an approximate 
synthesis of a control law containing delta 
functions is obtained by multiplying u§ by a 
large positive constant K. The total control 
is then given by the sum of the two components, 
with K appearing as a parameter in the 
feedback matrix. 


u 


-O.O86K 1.465-2K 2.875-1. 89K K 

O.O86K 1.465+2K 2.875+1.89K -K 


Figures 4(a) and 4(b) demonstrate the con- 
vergence properties of the control law as a 
function of the gain constant K. 


Although optimal control theory offers the 
most general method available for the design of 
model following systems, it is inefficient, 
because of the computational effort required, 
for answering such preliminary design questions 
as whether or not it is possible to match model 
and plant and whether bounded or unbounded con- 
trols are required. The theory presented here 
answers such questions directly by means of an 
algebraic test and, in addition, furnishes, for 
the case of implicit model following, a simply 
computed control law that achieves perfect 
matching if the test shows this to be possible. 
By means of an example it is demonstrated that 
even if perfect matching is not possible the 
performance of a system using the simple con- 
trol law may compare favorably with the perfor- 
mance of systems designed via optimal control 
as long as the dynamics of plant and model are 
not too dissimilar. 


Appendix 


expression for the time response of a linear 
system [5], we compute y(At): 


Approximate Synthesis of Feedback Lav 
Containing Delta Functions 

In the derivation of conditions for per- 
fect following it was necessary to include 
delta functions as permissible controls in 
order to achieve perfect following. Whenever 
the model is of lower order than the plant , it 
may be expeditious to sacrifice a part of the 
plant dynamics for better matching, and in that 
case delta function controls are necessary. 
Important questions arise now as to the proce- 
dure for constructing a control law containing 
delta functions and how to approximate one to 
arbitrary accuracy. Although the concept 
behind the construction of such a control law 
is well known, its adaptation to this problem 
requires some explanation. 


To begin with, it is assumed that perfect 
matching in the implicit model following sense 
can be achieved with controls containing finite, 
delta function and derivative of delta function 
components. These three components of u are 
given in the main text and are repeated here 
for convenience: 

Ui = Piu , U 5 = P5U , u^ = P gl u 

(Al) 

where u is defined as follows 


u = M T (LH - KF)x 


= (HB + HFBPg + HF 2 BP &1 )(LH - HF)x (A2) 


The next step is to divide time into equal 
increments At, which are chosen much shorter 
than the shortest time constant of the model 
and the plant- The control applied to the 
plant remains constant throughout each time 
increment and is updated only at the beginning 
of a new increment. Assuming the control pro- 
cess starts at t = 0, the first control 
applied to the plant is chosen as follows: 


u(o) = 


+ 5t p s 


(At) 2 V 




(LH - HF)x(o) 
(A3) 


We note that the gain constants multiplying 
those components of control that require delta 
function and derivative of delta function are 
l/At and 2/(At) 2 , respectively. Also, x(o) is 
arbitrary and therefore may be such that 
y(o) i Ly(o) . 


y(At) = H e FAt x(o) + 


Ft 


(AM 


where the transition matrix e is given in 
terms of the infinite series as follows: 


Ft _ , , F 2 t 2 , 

e — I + Ft + "2 1 + 


(A5) 


Upon substitution of equation (A5) into (kb) 
and since u(t) is constant within the integra- 
tion interval, y(At) becomes 

y(At) = H -jx(o) + AtFx(o) + ~ y F 2 x(o) + . . . 


IAt + 


FAt" 


F 2 At 3 

3-’ 


.Jbu(o' 


(a6) 


The purpose of calculating y(At) is to evalu- 
ate the error between y(At) and Ly(At) and to 
show that it can be made arbitrarily small. 

The error, denoted by ri(At), is evaluated 
using equations (l), (2), (A3), and (A6): 


7}(At) 

= H /f + AtF 2 + 


At 2 F 3 


FAt + 


F 3 At 3 


2BP. 


f~ + • • • + [ J 
r - • • -][ Bp i + It 

:] M t (EH - HF) + b|pi + i P S 
^Jm + (LH - HF)j- x(o) 


F 2 At ; 


(At) 2 
2 

(At) 2 ^ 


{ L 


- -{LH + AtLHP + ^|y LHF 2 + . . 


+ LH 


IAt + ^ + ?" At3 


r][ BPl 


BPs 

At ' (At) 


(At) 2 J 


(LH - HF ) V x(o) 


(A8) 


Since HBPg = 0, HBP 5 i = 0, and HFBP &1 = 0, 
equation (A8) simplifies to the following form: 

f|(At ) = (MM f - I)(LH - HF) + o(At) (A9) 


This sets the stage for the crucial step 
of this approach, namely the computation of the 
error between y and Ly at the end of the 
first time increment. Using the standard 


where o(At) are terms that go to zero at least 
as fast as At. But the first term is identi- 
cally zero by the assumption that the perfect 
matching condition, equation (l6), is satisfied. 


Thus, if At is chosen sufficiently small (or 
the gain constants arbitrarily large) , the 
error at the end of At seconds can he made as 
small as desired for any initial condition. 

Furthermore , the error can he maintained 
arbitrarily small for all future time if the 
control, equation (A 3 ), is updated at the 
beginning of each new time increment by replac- 
ing x(nAt) with x[(n + l)At] , n being the 
number of time increments. 

Through appropriate limiting arguments, 
the discussion given here for a discrete time 
control law can be suitably generalized to 
continuous time control. 
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Figure 1.- Transient responses pertaining to 
example 1 : initial bank angle 

disturbances. 



TIME, sec 

Figure 3- _ The effect of initial bank 
misalinement between model and plant 
in real model following. 


Figure 4 . - The effect of K on the trans ient 
response of example 2. 



